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Abstract 

By using a renormalization group analysis, we study the effect of interpar- 
ticle interactions on the critical temperature T^kt at which the Berezinskii- 
Kosterlitz-Thouless (BKT) transition occurs for Bose-Einstein condensates 
loaded at finite temperature in a 2D optical lattice. We find that Tbkt 
decreases as the interaction energy decreases; when U/J = SG/vr one has 
Tbkt = 0, signaling the possibility of a quantum phase transition of BKT 
type. 

I. INTRODUCTION 

It has been recently suggested [1] that a 2D optical lattice of Bose-Einstein condensates 
[2] may allow for the observation of a finite-temperature phase transition to a superfluid 
regime where the phases of the single- well condensates are coherently aligned. In fact, in 



an appropriate range of parameters, the thermodynamical properties of the bosonic lattice 
at finite temperature may be well described by the Hamiltonian of the XY model [1], and, 
as it is well known, the 2D XY model exhibits the Berezinskii-Kosterlitz-Thouless (BKT) 
transition [3-5] . The BKT phase transition occurs via an unbinding of vortex defects: the 
low-temperature phase, for T below the BKT temperature Tbkt, is characterized by the 
presence of bound vortex-antivortex pairs and the spatial correlations exhibit a power-law 
decay. For T ~ Tbkt the pairs starts to unbind (see e.g. [6,7]): in the high-temperature 
phase, only free vortices are present, leading to an effective randomization of the phases and 
to an exponential decay of the correlation functions. 

The Hamiltonian describing the properties of bosons in deep optical lattices is the so 
called Bose- Hubbard Hamiltonian [8], in which two terms are present: a kinetic term de- 
scribing the hopping of the bosons with tunneling rate t, and a potential term describing 
the interaction between bosons in the wells of the array with interaction energy U (which 
is proportional to the s-wave scattering length) . The description of bosons in deep lattices 
by means of the Bose-Hubbard Hamiltonian holds also at finite temperature, provided that 
the temperature effects do not induce the occupation of higher bands: for a 2D optical lat- 
tice confined to a plane by a magnetic potential having frequency a;^, this implies that the 
Bose-Hubbard Hamiltonian description is valid up temperatures T such that fiuiz ^ ksT. 
When the average number of bosons A^o per site is high enough and the interaction energy 
is larger than t/NQ, one may map - also at finite temperature - the Bose-Hubbard Hamil- 
tonian in the quantum phase Hamiltonian [1]: the conjugate variables are the phases and 
the particle numbers of the condensates in the different sites of the lattice. Thus, also at 
finite temperature, the phase diagram is determined by the competition of two energies: the 
Josephson energy J ^ 2tiVo, proportional to the tunneling rate between neighbouring sites 
of the 2D square optical lattice, and the interaction energy U. For U — 0, the quantum 
phase Hamiltonian reduces to the XY Hamiltonian, exhibiting then a BKT transition at 
Tbkt ~ J/ks'- for T < Tbkt the system as a whole behaves as a supcrfiuid with phase coher- 
ence across the array. Very accurate Monte Carlo simulations yield, in the thermodynamic 
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limit, Tbkt = 0.898 J/Zcb [9]. 

In this paper we determine the effect of the interaction energy U on the BKT transition. 
When U ^ J, Si BKT transition still occurs at a critical temperature Tbkt{U): we find 
that Tbkt{U) is smaller than Tbkt and decreases with U (see Fig. 1). Intuitively speaking, 
the quantum phase Hamiltonian, as well as the Bose-Hubbard Hamiltonian, in the limit 
U ^ J describes a Mott insulator, while in the opposite limit, J ^ U, the array behaves 
as a phase-coherent superfluid (e.g., see [10,11]): thus, when U increases, the superfluid 
region in the phase diagram decreases. Although in two dimensions there is not long-range 
order at finite temperature due to the Mermin- Wagner theorem, still the system exhibits 
superfluid behavior for T < Tbkt [H]: thus one expects that Tbkt{U) should decrease 
when U increases. We also find (see Section III) that, when [//J is equal to the critical 
value {U/J)cr = 36/7r, one has that Tbkt{U) = 0, signaling the possibility of a (T = 0) 
quantum phase transition. 

The quantum phase Hamiltonian describes also the behavior of superconducting Joseph- 
son networks below the temperature Tbcs at which each junction becomes superconducting 
[11-13]. The momenta conjugate to the phases of the macroscopic wavefunctions of the 
superconducting grains are the number of Cooper pairs, J is the Josephson energy of the 
superconducting Josephson junctions and U is the charging energy due to the Coulomb re- 
pulsion between Cooper pairs. We observe that, in superconducting Josephson arrays, the 
interaction term is written in general as J2ij UijNiNj where Ni is the Cooper pairs number 
at the site i: Uij is proportional to the inverse of the capacitance matrix and may be also 
non-diagonal (i.e., Uij ^ with i ^ j). At variance, for bosons in deep optical lattices, 
the interaction term is J2i UuN^, which corresponds, in a suitable range of parameters, to a 
diagonal quantum phase model. 

The analogy between superconducting Josephson networks and atomic gas in deep optical 
lattices is then clear: in each well of the periodic potential there is a condensate grain, 
appearing at the Bose-Einstein condensation temperature Tbec- When Tbec is larger than 
all other energy scales, the atoms in the well i of the 2D optical lattice may be described 
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by a macroscopic wavefunction. It becomes apparent, then, that an optical network can be 
regarded as a network of bosonic Josephson junctions. Furthermore, in both systems, when 
the interaction term is neglected respect to the energy associated to the particle hopping, 
one has that the whole 2D array becomes superfluid at the temperature Tbkt at which 
the BKT transition occurs, with Tbkt smaller than Tbcs for superconducting Josephson 
networks and smaller than Tbec for bosonic Josephson networks. 

The plan of the paper is the following: in Section II we introduce the effective Hamil- 
tonian describing bosons hopping on a deep optical lattice and, by using a semiclassical 
approximation [14], we compute the effective Josephson energy in presence of the interac- 
tion energy U. In Section III we evaluate the effect of the quantum fluctuations on Tbkt by 
putting in the renormalization group equations the effective Josephson energy obtained in 
Section II; a comparison with previous results is then carried out. Details of the computa- 
tions skipped in the text are in the Appendix A. Section IV is devoted to some concluding 
remarks. 



2D optical lattices are created using two standing waves [2]: when the polarization 
vectors of the two laser fields are orthogonal, the periodic potential is 



where k — 2tt/X is the wavevector of the lattice beams. Vq is usually expressed in units of 
Er — Tb^k'^ /2m (where m is the atomic mass): in [2] it is A = 852 nm and Er — h ■ S.lAkHz. 
Around the minima of the potential (1) one has Vopt ~ muj^{x'^ + y^)/2 with 



When ujr ^ ujz (with cu^ the frequency of the confining magnetic potential superimposed to 
the optical potential and acting along z), the system realizes a square array of tubes, i.e. an 
array of harmonic traps elongated along the 2;-axis [2] . 



II. THE RENORMALIZED JOSEPHSON ENERGY 



Vopt = Vb[sin2 (kx) + sin^ (ky)], 



(1) 




(2) 
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When all the relevant energy scales are small compared to the excitation energies, one 
can expand the field operator [8] as 

if{r,t)=J2^,{t)Ur) (3) 

i 

with $i(r) the Wannier wavefunction localized in the i-th well (normalized to 1) and Ni = 
'ipl'ipi the bosonic number operator. Substituting the expansion of ^(r, t) in the full quantum 
Hamiltonian, one gets the effective Hamiltonian describing the bosons hopping on the deep 
optical lattice [8,15] 

H^-tY. (^14- + h.c.) + ^Y.HN^ - 1). (4) 

<i,j> i 

In Eq.(4) Y.<i,j> denotes a sum over all the distinct pairs of nearest neighbours; t and U are 
respectively the tunneling rate and the interaction energy, and are given by 
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t C::i — j df 

and 



h 

2m ^ I ext 3 



(5) 



U 



m 



J dr^l (6) 



(a is the s-wave scattering length). 

As discussed in the Introduction, one can show that, when Vq and ujz are large enough, 
the system is described by the Hamiltonian (4) up to temperatures T ~ hcUz/kB [!]• Upon 
defining J pa 2tNo, the Hamiltonian (4), for Nq 1 and J/Nq <^ U [16], reduces to 

H=-J E cos(^,-^,)-|x:^. (7) 
<i,j> 2 . 

For a 2D lattice with Vq between 20Eji and 25Eii and Nq pa 170 as in [2], one sees that the 

condition » J/Nq is rather well satisfied and that J/ks is of order of 20nK. For U — 0, 

(7) is the XY Hamiltonian: 

H^-J cosiOi-Oj). (8) 
<hj> 



The Hamiltonian (7) describes the so-called quantum phase model [11-13]. There is an 
huge amount of literature on the properties of 2D superconducting Josephson arrays studied 
by means of the quantum phase model [11-13]. A (not exhaustive) list of relevant papers in- 
cludes mean-field and coarse-graining approaches [17-19], Monte Carlo results [20-24], renor- 
malization group calculations [14,21] and self-consistent harmonic approximations [21,25,26] 
(more references are in [11]). In the following we shall study the renormalizat ion-group 
equations in which it is used an effective value of the Josephson energy computed within the 
harmonic approximation [14]. 

The starting point is the partition function Z of the quantum phase model: using the 
path integral formalism, from (7), one has 



z ^ j vee-"n'^'\ 



(9) 



where the action S is given by 

H 

s[e] ^ JdT 



J \ / <2J> 



(10) 



with P — l/ksT and 9ij = 9i — 9j. Separating the phases as ^^(r) — (pi + ^^(t), where ipi is 
a static vortex configuration and ^^(t) is a fluctuation about (fi, the path-integral partition 
function (9) can be written as 



(11) 



where 



S[ip, S] = J 



rt.2 



dr 



+ J (1 - cos (fiij COS Sij) 

i <i,j> 



(12) 



with (fiij = (fi — ipj e dij = 5i — Sj. Assuming that ipi and 5i are slowly varying over the size 
of the array [14,12], i.e. 



cos ipij 



(13) 



and 



COS 5ij Ki 1- 5'fJ2, 



(14) 



one gets 



where Zq = JVSe-^^°^^^ and 



<i,j> 



hp 

So[5] = Jdr 





■ J 

y 52 ^ :^ y ^2 

I <i,j> 



In Eq.(15) J is the renormahzed Josephson energy, which is given by 



j(i-^<4>o) 



with 



The evaluation of < 5fj >o can be carried out in a standard way [12] and one has 



< >o= Y -J- —^jx^ coth X dx 



where 



(15) 



(16) 



(17) 



(18) 



(19) 



ri = pV7rUJ. 



(20) 



III. RENORMALIZATION GROUP EQUATIONS 

In the renormalization group equations for the 2D XY model [3-5] the scale-dependent 
screened charge /C depends on the dimensionless scaling parameter I = log(r/a) (where r 
is the vortex distance and a is the lattice spacing) and it is given by /C(/) = pj/e{l), where 
the dielectric constant e{l) expresses the screening of the vortex- antivortex interaction due 
to the presence of other vortices [5] . The renormalization group recursion equations read [5] 
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dl 



(0 



47rV(/) 



(21) 



and 



dy{l) 



[2 - n)C{lMl) 



(22) 



dl 



where y oc r 



•2g /3y(r)/2g f3n^ with 6 is the fugacity for creating a vortex pair and V{r) 



correspond to the screened vortex-antivortex potential. The study of the scahng equations 
(21)-(22) about the fixed point y/ — 0, }Cf — 2/7r shows that the BKT transition occurs at 
2 - 7r}C(l = 0)^0, where }C{1 = 0) = (3J/e{l = 0)=(3J [5]. 

In presence of the quantum fluctuations, one has to replace the scale-dependent charge 
/C(0 by jC{l) with ^(Z = 0) = /3J and J given by Eqs.(17) and (19). Denoting with 
Tbkt{U) the BKT transition temperature for a given [/, one finds the following equation 
ior K=J/kBTBKT{U): 



where = U/ttJ. 

The root of Eq.(23) indicates the critical point at which a BKT transition occurs. The 
eventual occurrence of a double root for Eq.(23) might correspond to what is called in 
literature a reentrant behavior. Indeed it has been often argued (see Refs. in the review 
[11]) that the quantum phase model may undergo at low temperatures a reentrant transition 
induced by the quantum fluctuations: namely, fixing U / J and lowering the temperature, 
one could switch from an insulating phase to a superconducting one at Tbkt{U) and then, 
lowering further the temperature, one finds another critical temperature T^^\U) at which 
the system comes back {reenters) in the insulating phase. Consistent with the reentrant 
scenario is the dramatic decrease of the specific heat at very low temperatures [20] and, as 
discussed, the presence of double roots for Tbkt in the renormalization group equations. 

The phenomenon of the reentrance - although not universal - is a non-perturbative effect: 
in fact, opposite results may be obtained by means of different truncations for the series 




(23) 
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expansion of the function F. For instance, if in Eq.(23) one expands the hyperbohc cotangent 
as cothx 1/x + x/3, one gets an expression for F{K) which gives two roots for Tbkt- If 
instead, for r? < tt, one uses the expansion cotha; = - + Kl^?,"- x^'^~^ where the B^s are 

n=l ^ ' 

the BernouUi numbers [27], one obtains the equation: 



= 



(24) 



with Xu — 7rC//24J. Since B/^n < for n = 1, 2, . . ., and -B4„+2 > for n = 0, 1, 2, . . ., if 
one truncates the sum in Eq.(24) to the (2n + l)-th order with n — 0, 1, 2, . . ., one finds 
F{K) ^ oo for — > oo and then Eq.(24) has two roots; at variance, if one truncates the 
sum to the (2n)-th order with n = 1, 2, . . ., one finds F{K) — > — oo for X — > oo and Eq.(24) 
admits only one root. Thus one cannot truncate the sum in Eq.(24) to any finite order, even 
for U/J <C 1; the origin of this problem is that the expansion used to get (24) is a series 
with terms having alternating signs. In addition, to find the roots of the equation F{K) — 
one has to evaluate the integral in (19) up to \/UJ /ksTsKT, and, also with U small, 77 may 
become large requiring to use all the orders in Eq.(24). 

A way to overcome the above mentioned difficulties is to use in Eq.(23) the Mittag-Leffler 

00 

expansion of the hyperbolic cotangent, i.e. vrcothTrx = l/x-\-2x (x^ + n^)~^: in this way 
Eq.(23) can be written as 

1 1 



n=l 



F{K) = 2-TiK 



1 - 







(25) 



where 



g{K, X^) = X; - Y log(l + 



n=l 



(26) 



An analytic expression for F{K) valid for f// J ^ 1 is given by Eq.(A5) in the Appendix A. 
A detailed study of the function F{K) is in the Appendix A. Here we observe that, in 



the limit U/J 0, one finds 2 + f - ttX = 0, from which Kq = K{U/J ^ 0) = ^ and 
ksTEKTiS^) — 1.128 J. Furthermore, for U/J < 36/tt, one can show (see the Appendix A) 
that F{K) — > —00 for — > 00 and that F'{K) < 0: since F tends to the positive value 
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2 + 7r/4 for X — > 0, one can safely conclude that, for U/ J < 36/7r, Eq.(25) has an unique 
solution. At variance, one can show that, for U/J > 36/7r, Eq.(25) does not admit any 
solution, while for U/J = 36/7r one has K ^ oo, i.e. Tbkt{U) — > 0. A plot of F{K) 
and F'{K) for three values olU/ J respectively smaller than, equal to and larger than 36/7r 
is given in Figs. 2 and 3. One may infer that, at T = 0, a phase transition is expected 
to occur at the critical value {U/J)cr = 36/7r: this value is in reasonable agreement with 
the mean-field estimates for the T = transition. In fact, for the diagonal quantum phase 
model, one has that the mean-field prediction [17,19] is {U/J)cr ^ 2z — % where 2; = 4 is 
the number of nearest neighbours. 

In Fig. 1 we plot as empty circles the values of Tbkt{U) /Tbkt{^) as a function of C//J 
from the numerical solution of Eq.(25). If one uses the analj^ic expression (A5) for the 
function F{K) one gets, for t//J ^ 1/2, an error lesser than 1%. A much better estimate 
may be obtained by expanding the function F{K) near Kq: in this way Eq.(25) reads 

F{K) ^ F{Ko) + {K- Ko) ■ F\K^) = 0, 

and then 

^ ~ ^° " n^y ^^^^ 

In Fig. 1 Eq.(27) is plotted as a sohd line. 

We may conclude that Eq.(25) admits, for small C// J, a unique solution and this excludes 
the possibihty of reentrant behavior. Of course, our conclusion rehes on the approximations 
(13)-(14) made in order to estimate J and < Sf^ >q. A more careful treatment accounting for 
the periodicity of the phases in Eq.(12) seems to be highly desirable to put on a more solid 
base all the contrasting issues related to the reentrant behavior of the systems described 
by the Hamiltonians (4) and (7). It is comforting to see that the experimental study of 
weakly interacting bosons on 2D optical networks and the investigations of their superfluidity 
properties at very low temperatures could provide new insights also on this intriguing and 
yet poorly understood problem (see also the recent papers [28]). 
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IV. CONCLUDING REMARKS 



In this paper we studied, by means of a renormalization group analysis, the effect of 
interparticle interactions on the critical temperature at which the Berezinskii-Kosterlitz- 
Thouless transition occurs for Bose-Einstein condensates loaded in a 2D optical lattice at 
finite temperature. We determined the shift of the Berezinskii-Kosterhtz-Thouless transition 
temperature induced by the interaction term and we compared our findings with previously 
known results. 

APPENDIX A: PROPERTIES OF THE RENORMALIZATION GROUP 

EQUATION 

In this Appendix we study the properties of the function F{K) defined in Eq.(25). We 
show in (1) that F(0) > 0, in (2) that F{K) -oo for X ^ oo when U/J < 36/7r and 
in (3) that for C//J < 36/7r one has F'{K) < 0: one may then conclude that the equation 
F{K) — has only one root for U/J < 36 /tt. Finally, we derive an analytic expression for 
F(K) holding for small U/J. 

Prom Eq.(25) one can see that: 
(1) for U/J ^ one has 2 + f - tt/T = 0, from which K = ^ and KbTbkt ^ 1.128J: 
indeed, for small U/J, one has that 

2 2 ^ 2 2 V n2 2n^ ^"^ 

and thus 

In the same way, F{K) ^ 2 + f for X ^ 0. 

(2) F{K) ^ -oo per /s: ^ oo for Xu < 3/2 (where = 7r^Xj24 = 7rU/24J): indeed 
for large K 
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" 

then, for K — > oo, F[K) — > — oo with < 3/2 and F[K) — > oo with x^ > 3/2. 



(3) For Xu < 3/2, it is F'{K) < for /sT > 0: indeed 

F'{K) = -TT - g{K, X^) + [coth(7rXy^i^) - — ^] . (Al) 



Since hmj-^o (coth a; ~ ^) = 0, one gets 



hm F'iK) = -TT. 



For large K one finds 



Therefore one has that \\m.K^ooF' [K) < for < SG/tt"^, and hm/^^oo > for 
Xu > 36/7r^. Using similar arguments, one can show that, for K > 0, it is 

F'{K)<-n + ^Xl/\ (A2) 

which implies that, for X„ < 36/7r2, one has F\K) < 0. The behavior of F{K) and F\K) 
for three values oiU/J respectively smaller than, equal to and larger than 36 /n is plotted 
in Figs. 2 and 3. 

One may also obtain an analytic approximation for F{K) holding for [// J ^ 1 by putting 
z — \fXuK and using log(l ^ z) — ■ ^ for | z |< 1: one has 



5(i^,X„) = EW— C(2i-2) (A3) 

j=2 



oo 

where C,{j) = ^ is the Riemann zeta- function. One finds 

n=l 
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For large n one has the recurrence relation [27] 



C(n + l)^-[l + CW]. 



Substituting back in Eq.(A3) and using (A4), one has 



g{K, X^) ^ X^K' (y - 2) + 2 log J^.^, ^ + ^ log(l + X^K') 



siiihTiXi' K 2 



and finally 
F fa 2 - ttX 



(A5) 
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FIGURES 




FIG. 1. The BKT critical temperature Tbkt{U) (in units of Tbkt{'^)) as a function of U/J. 
Empty circles: numerical solution of Eq.(25); solid line: Eq.(27). 




FIG. 2. F{K) for U/J = 1 (solid line), 36/ tt (dashed line) and 15 (dotted line). The dot-dashed 
line is the analytic approximation (A5) for F holding for U/J<^ 1. For C//J = 36/7rF asymptot- 
ically tends to 0. 



16 



- . ,.- 



■2 j 



5 ^ 7« 



FIG. 3. F'{K) for [//J = 1 (solid line), 36/7r (dashed line) and 15 (dotted line). 
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